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Abstract 

We describe in detail an improved standard solar model which 
has been used to calculate the fluxes of standard solar neutrinos. It 
includes premain sequence evolution, element diffusion, partial ioniza- 
tion effects, and all the possible nuclear reactions between the main 
elements. It uses updated values for the initial solar element abun- 
dances, the solar age, the solar luminosity, the nuclear reaction rates 
and the radiative opacities. Neither nuclear equilibrium, nor com- 
plete ionization are assumed. The calculated solar neutrino fluxes are 
compared with published results from the four solar neutrino experi- 
ments. The calculated solar neutrino flux is consistent, within the 
theoretical and experimental uncertainties, with the solar neutrino ob- 
servations at Homestake and Kamiokande. The observations suggest 
that the '^Be solar neutrino flux is much smaller than that predicted. 
However, conclusive evidence for the suppression of the ^Be solar neu- 
trino flux will require experiments like BOREXINO and HELLAZ. If 
the ''Be solar neutrino flux is suppressed, it still can be due either to 
standard physics and astrophysics or neutrino properties beyond the 
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standard clcctroweak model. Only future neutrino experiments, such 
as SNO, Superkamiokande, BOREXINO and HELLAZ, will be able 
to show that the solar neutrino problem is a consequence of neutrino 
properties beyond the standard electroweak model. 
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1 Introduction 



The Sun is a typical main sequence star which generates its energy 
by fusion of protons into Hehum nuclei through the pp and CNO 
nuclear reactions chains (Bethe 1939). These reactions also produce 
neutrinos. These neutrinos have been detected on Earth in four pio- 
neering solar neutrino experiments, the radiochemical Chlorine exper- 
iment at Homestake (Cleveland et al. 1995 and references therein), 
the electronic light water Cerenkov experiment at Kamioka (Suzuki et 
al. 1995 and references therein) and the two radiochemical Gallium 
experiments, GALLEX at Gran Sasso (Anselmann et al. 1995 and ref- 
erences therein) and SAGE at the Baksan (Abdurashitov et al. 1995 
and references therein). They provide the most direct evidence that 
the sun generates its energy via fusion of Hydrogen into Helium. How- 
ever, it has been claimed (e.g., Bahcall 1995, Hata et al. 1995) that all 
four experiments measured solar neutrino fluxes significantly smaller 
than those predicted by standard solar models (SSM) (e.g., Bahcall 
and Uhich 1988; Turck Chieze et al. 1988; Sackman et al. 1990, 
Bahcall and Pinsonneault 1992 (BP92), Turck-Chieze and Lopes 1993 
(TL93); Castellani et al. 1994; Kovetz & Shaviv 1994; Christcnsen 
- Dalsgaard 1994; Shi et al. 1994; Bahcall and Pinsonneault 1995 
(BP95); See, however, Dar and Shaviv 1994 (DS94); Shaviv 1995; 
Dzitko et al. 1995. 

The Homestake chlorine experiment has been collecting data 
since 1970. Based on 25 years of observations they have recently 
reported (Cleveland et al. 1995) an average ^^Ar production rate 
in 3^01 of 

Rate (Chlorine) = 2.55 ± 0.17(stat) ± 0.18(syst) SNU, (1) 

(where {ISNU = 10"^^ captures per atom) by solar neutrinos 
above the 0.814 MeV threshold energy for the reaction i/g-l-^^Cl— > 
e-+37Ar. It is 32 ± 5% of the 8.lt\i SNU predicted by the SSM of, 
e.g., Bahcall and Pinsonneault 1992 (BP92). 

Kamiokande II and III observed electron recoils, with energies 
first above 9 MeV and later above 7 MeV, from elastic scattering of 
solar neutrinos on electrons in water. Their 5.4-year data show a 
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spectrum consistent with solar neutrino flux of (Suzuki 1995) 



= [2.9 ± 0.2{stat) ± 0.3{syst)] x 10^ cm-^g-i , (2) 

which is 51% ± 9% of that predicted by the SSM of BP92. 

GALLEX, the European Galhum experiment at the Gran Sasso 
underground laboratory, measured a ''^Gc production rate by solar 
neutrinos in ^^Ga of (Anselmann et al. 1995a) 



through the reaction fe+ Ga— > Ge whose threshold energy is 

233 KeV. This rate is 60% ± 10% of the 131.5^^ SNU predicted by 
the SSM of BP92. 

SAGE, the Soviet-American Gallium Experiment in the Baksan 
underground laboratory reported (Abdurashitov et al. 1995) an aver- 
age ^^Ge production rate by solar neutrinos in ^^Ga of 



during 1990-1993, which is consistent with the rate reported by GALLEX. 
It is 56% ± 11% of the 131.51^ SNU predicted by the SSM of BP92. 

These discrepancies between the observations and the predictions 
have become known as the solar neutrino problem(s). Three types of 
solutions to these solar neutrino problems have been investigated: 

(a) Terrestrial Solutions: Perhaps the accuracy of the results of 
the solar neutrino experiments has been overestimated and unknown 
systematic errors are largely responsible for the solar neutrino prob- 
lem. 

(b) Astrophysical Solutions: Perhaps the standard solar models 
do not provide sufficiently accurate description of the present day sun 
and/or the neutrino producing reactions in the sun. 

(c) Particle Physics Solutions: Perhaps non standard neutrino 
properties beyond the standard model are responsible for the solar 
neutrino problem. 

The chances that possibility (a) is responsible for the solar neu- 
trino problem have been greatly reduced by the GALLEX Chromium 
source experiment performed last year (Anselmann et al. 1995b). 
This experiment is the first full demonstration of the reliability of the 



Rate (Gallium) = 79 ± lO(stat) ± 6(syst) SNU, 



(3) 



Rate (Gallium) = 74tllistat)tr{sys) SNU 



(4) 
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radiochemical technique for the detection of solar neutrinos. In partic- 
ular, it excludes the possibility of any unidentified important sources 
of systematical errors, such as hot atom chemistry, in the radiochem- 
ical experiments. However, we will argue that terrestrial solutions 
to the solar neutrino problem have not been ruled out. In particu- 
lar, the capture rate of solar neutrinos near threshold in ^^Cl and in 
^^Ga may be significantly smaller than currently used for predicting 
SNUs (see chapter 7). We also note that although Kamiokandc, the 
electronic water Cerenkov experiment, has definitely observed excess 
events pointing towards the sun (Suzuki 1995), it has not been cali- 
brated yet with a known neutrino source. Such a direct calibration is 
highly desirable since Kamiokande has the formidable task to single 
out solar neutrino events with an expected rate of the order of 1/day 
in a fiducial volume of 680 tons of highly purified water where the 
total trigger rate of the detector for electron energy greater than 7 
MeV is still 150,000/day 

It should also be noted that the joint observations of Homestake 
and Kamiokande appear to be supported by the results of the Gallium 
experiments (Dar 1993). Namely, both the joint results from Homes- 
take and Kamiokande and the results from GALLEX and SAGE seem 
to indicate that the '''Be solar neutrino flux is much smaller than that 
predicted by the SSM (e.g., Dar 1993; Hata et al. 1994; Berezin- 
sky 1994; Kwong and Rosen 1994; Parke 1995). In fact, the signals 
observed by GALLEX and SAGE can be predicted accurately from 
the observed solar luminosity and the solar neutrino observations of 
Homestake and Kamiokande (Dar 1993): 

If the sun derives its energy from fusion of Hydrogen into Helium 
and if it is in a steady state where its nuclear energy production rate 
equals its luminosity, then the total solar neutrino flux at Earth is 
given by (e.g., Dar and Nussinov 1992) 



6.51 X 10^° cm~2g-i ^ 



where D 1.496x 10^^ cm is the distance to the sun, Q = 26.733 MeV 

is the energy released when four protons fuse into Helium, and the 
average energy of solar neutrinos has been approximated by Eiy{pp) = 
J2 Ey^i(j)y^i (piy^i « 0.265 M cV , the average energy of the pp neutrinos 
which dominate the sum rule. If the pep and CNO contributions to 
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the solar neutrino flux are ignored then, 



{pp) + ^.^{^Be) + <l),^{^B). 



(6) 



Kamiokande reported 



^^^(85) = (2.9 it 0.4) X 10^ cm-^s" 
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(7) 



Thus, the flux of B neutrinos alone contributes < cr^ >= 3.2 it 
0.50 SNU to the '^'^ Ar production rate in ^'''Cl, assuming a cross sec- 
tion of (1.11 ± 0.03) X 10~^^cm? for the capture of '^B neutrinos by 
^^Cl (Aufderheide et al. 1994). The Chlorine experiment at Homes- 
take which is sensitive also to ^Be neutrinos observed since 1970 an 
37 Ar production rate of only 2.55 ± 0.25 SNU, and 2.78 ± 0.35 SNU 
in runs 91-124 during 1986-1993, the running period of Kamiokande 
II and III. For standard solar neutrinos these results leave very little 
room for a ''Be solar neutrino flux (the SSM estimates its expected con- 
tribution at a level of ~ 1 SNU). Therefore, we conclude (Dar 1993) 
that Be) ^ Be) and (pjy^ « (pu^ipp)- Consequently, the 

expected ^^Ge production rate in ^^Ga is given by 



< <T(/>^0 >Ga~< (^(I>uq{Lq) >Ga + < acji^^iKAM) >Ga= 83 SNU, 



where the first term on the rhs is the contribution of the pp solar 
neutrino flux estimated from the solar luminosity and the second 
term is the contribution from the ^B solar neutrino flux measured by 
Kamiokande. This predicted signal is in excellent agreement with the 
79 ± 12 SNU measured by GALLEX and the 73 ± 19 SNU measured 
by SAGE, respectively. 

Note, however, that in view of the experimental uncertainties, the 
observed signal in the Chlorine experiment is not signiflcantly smaller 
than the expected minimal signal due to the solar neutrino flux mea- 
sured by Kamiokande. Thus, the solar neutrino experiments do not 
provide conclusive evidence for new electroweak physics as claimed, 
for instance, by Bahcall and Bethe (1991). Although neutrino oscilla- 
tions, and in particular the MSW effect (Mikheyev and Smirnov 1986, 
Wolfenstein 1978) can solve the solar neutrino problems, the inferred 
neutrino mixing parameters seem to differ substantially from those im- 
plied by the atmospheric neutrino anomaly (for a recent summary sec 
e.g., Barish 1995) or those implied by the neutrino anomaly observed 



(8) 
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by the LSND collaboration at LAMPF (Louis 1995; Atlianassopolous 
et al. 1995; see, however Hill 1995) . Moreover, in spite of extensive 
laboratory searches and many precision tests no confirmed evidence 
from accelerator experiments has been found for new physics beyond 
the standard electroweak model (e.g., Langacker 1995 and references 
therein). Without such an evidence it is quite important to improve 
the standard solar model (which provides only an approximate de- 
scription of the complicated sun) and to continue the search for a 
standard physics solution to the solar neutrino problems. 

In this paper we present an improved standard solar model and its 
predictions for the solar neutrino fluxes. The model includes premain 
sequence evolution, diffusion of all elements, partial ionization effects, 
and all the significant nuclear reactions between the various elements 
which the sun is made of. It uses updated values for the initial solar 
element abundances, the solar age, the solar luminosity, the nuclear 
reaction rates and the radiative opacities. Neither nuclear equilibrium, 
nor complete ionization are assumed. It employs a very fine zoning of 
the sun and accurate numerical procedures to integrate the solar evo- 
lution equations from zero age until the present day. The calculated 
solar neutrino fluxes are compared with those measured by the four 
solar neutrino experiments. The calculated solar neutrino flux is 
consistent, within the theoretical and experimental uncertainties, with 
the solar neutrino observations at Homestake and Kamiokande. How- 
ever, the observations appear to suggest that the ^Be solar neutrino 
flux is much smaller than that predicted by our solar model. In spite 
of the facts that the value of the ^Be solar neutrino flux is a robust 
prediction of the current standard solar models and that minimal ex- 
tensions of the standard electroweak model, such as neutrino flavor 
mixing, can explain its suppression, as will be explained in chapter 7, 
we do not consider the ^Be deficit to be a compelling evidence for new 
electroweak physics, as claimed by various authors (e.g., Bahcall and 
Bethe 1993; Bludman et al. 1993; Castellani et al. 1994; Hata et al. 
1994; Bcrczinsky 1994; Kwong and Rosen 1994; Bahcall 1994; Parke 
1995; Hata and Langacker 1995). This is because neutrino absorption 
cross sections near threshold in ^'''Cl and in '''^Ga may be significantly 
smaller than those calculated by Bahcall (1989), because the standard 
solar models arc only approximate and simplified descriptions of the 
real and complex sun, and because very little is known, either experi- 
mentally or theoretically, on dense plasma effects on nuclear reaction 
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rates, decay rates, particle and energy transport at solar conditions. 

Our paper is organized as follows: In Section 2 we outline the 
standard solar model. In Section 3 we describe the astrophysical pa- 
rameters which we use in this work. The main physics input is de- 
scribed in Section 4. Section 5 outlines our stellar evolution code. 
Our main results are described in Section 6 and are compared there 
with experimental results and the results of other solar model cal- 
culations. Possible standard physics solutions to the solar neutrino 
problem are briefly outlined in section 7. The neutrino mixing pa- 
rameters suggested by the MSW solution are derived analytically in 
section 8. Final conclusions are drawn in section 9. 

2 The Standard Solar Model - An 
Overview 

The standard solar model (e.g. Bahcall 1989 and references therein) 
is a physical description of the sun based on the standard stellar 
evolution equations, (e.g., Clayton 1968) which arc used to cal- 
culate its evolution from the premain sequence Hayashi phase to its 
present state, using the best available input physics (initial con- 
ditions, equations of state, nuclear cross sections, radiative opacities, 
condensed matter effects). The model assumes a complete spherical 
symmetry, no mass loss or mass accretion, no angular momentum gain 
or loss, no differential rotation and a zero magnetic field through the 
entire solar evolution . Thus, the assumed initial conditions are: 

1. Fully convective, homogeneous, spherically symmetric protostar. 

2. Initial mass of Mq = 1.99 x 10^^ gm. 

3. No angular momentum, no differential rotation, no magnetic 
field. 

4. Initial chemical composition deduced from primitive meteorites, 
the solar photosphere, the solar wind, the local interstellar gas 
and the photospheres of nearby stars. 

The calculations are iterated, treating the unknown initial ^He 

abundance and the mixing length in the convective zone (roTighly the 
size of the pressure scale height) as adjustable parameters, until the 
present day properties of the sun are reproduced. These include: 
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1. The solar luminosity Lq = 3.844 x lO^'^erg ■ sec ^. 

2. The solar radius Rq = 6.9599 x W^^cm. 

3. The observed solar surface element abundances. 

4. Internal structure consistent with helioseismology data (optional). 

The output of the calculations includes the present-day density 
profile p{r), temperature profile T{r) and chemical composition profile 
[Xj(r)] of the sun. They can be compared with information extracted 
from helioseismology. They can also be used to calculate the expected 
fluxes of solar neutrinos. In particular, according to the standard 
model, solar neutrinos are produced mainly in the fusion of hydrogen 
into deuterium {p + p ^ D + + Vg and p + e + p ^ D + u^), in 
electron capture by ^Be and in /? decay of ^B, ^^N, ^^O and ^^F. Their 
production rates in the sun are calculated using the standard elec- 
troweak theory and the density, temperature and element abundances 
in the sun provided by the standard solar model. Our calculations 

were performed with an updated version of the the solar evolution 
code of Kovctz and Shaviv. We refer the reader to its detailed de- 
scription by Kovetz and Shaviv (1994) and focus here only on main 
points and important improvements. 

3 The Astrophysical Parameters 
3.1 Solar Luminosity 

The absolute luminosity of the sun changes with time due to solar 
evolution. In addition to this secular solar variability the solar lumi- 
nosity changes with solar activity. The solar activity, its period and 
erratic disappearance (last time during the Maunder minimum) are 
not well understood. Precise measurements of the solar "constant" 
during a large fraction of an entire solar cycle have been carried out 
on a variety of spacecrafts and satellites during the last two solar cycles 
(21 and 22). They include Nimbus 7 spacecraft (see, e.g., Hickey et 
al. 1982), the Solar Maximum Mission Satellite (see, e.g., Wilson and 
Hudson 1988), The Earth Radiation Budget Satellite, the N0AA9 and 
NOAAIO satellites (see e.g., Lee et al. 1991) and The Upper Atmo- 
sphere Research Satellite (see e.g., Wilson 1993). Using the following 
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parametrization of the variation of the solar constant (SC) during the 
solar cycle (Wilson and Hudson 1988), 

SC{t) = S* X (1 + b{cos[2Tr{t - 1980.82)/10.95])) W ■ m'^ (9) 

we found a weighted mean average S = 1367.2 ± 3.4 W ■ for the 
solar constant. It yields a mean solar luminosity of Lq = AttD'^S = 
3.844 X 10^^ erg ■ s~^. This value is consistent with the value quoted 
by the Particle Data Group (1994) and used in BP95. It is smaller 
by 0.4% than the value Lq = 3.86 x 10^^ erg ■ s'^ used in BP92 and 
larger by 0.5% than Lq = 3.826 x 10^3 erg ■ s quoted before by the 
Particle Data Group (1992) and used by Dar and Shaviv (1994). 

3.2 Solar Age 

The solar age is best estimated from radioactive dating of mctcoritic 
condensation based on the observed 235^/238^ ^^^^ 207p)-,y206p)-, 
tios. Tilton (1988) estimated this age to be 4.56 x lO^y. Guen- 
ther (1989) estimated it to be 4.49 x lO^y while recent analyses by 
Gopel et al. (1994) and by Wasserburg et al. (1995) yielded Iq = 
(4.566 zb 0.005) X lO^y. This value was used by us as the total solar 
age which includes both the premain sequence evolution (some 30-40 
million years) and the time spent by the sun on the main sequence. 

3.3 Initial Chemical Composition 

The initial element abundances influence significantly the solar evolu- 
tion and the present density, chemical composition and temperature in 
the solar core, which determine the solar neutrino fluxes. In particular, 
the calculated radiative opacities, which in turn determine the tem- 
perature gradient in the solar interior, are very sensitive to the heavy 
elements abundances (the heavy elements are not completely ionized 
in the sun). There are four major sources of information on the initial 

solar abundances. They are: the chemical composition of the most 
primitive class of meteorites (type I carbonaceous chondrites), the so- 
lar photospheric abundances, the chemical composition of the solar 
wind and the local interstellar medium element abundances. 
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3.3.1 The Heavy Metals Abundances 



Apart from the noble gases, only a few elements, namely, H, C, N and 
O, have escaped complete condensation in primitive early solar system 
meteorites because they were able to form highly volatile molecules 
or compounds (see, e.g., Sturenburg and Holweger 1990). Thus, the 
initial solar abundances of all other elements arc expected to be ap- 
proximately equal to those found in type I carbonaceous chondrites 
as a result of their complete condensation in the early solar system. 
Since the chemical composition of the solar surface is believed to have 
changed only slightly during the solar evolution (by nuclear reactions 
during the Hayashi phase, by diffusion and turbulent mixing in the 
convective layer during the main sequence evolution, and by cosmic 
ray interactions at the solar surface) it has been expected that the 
photospheric abundances of these elements are approximately equal 
to those found in CI chondrites. Over the past decades there have 
been many initial disagreements between the meteoritic and photo- 
spheric abundances. In nearly all cases, when the atomic data were 
steadily improved and the more precise measurements were made, the 
photospheric values approached the meteoritic values. The photo- 
spheric abundances are now as a rule in very good agreement with the 
meteoritic values (Grevesse and Noels 1993a). Hence, as meteoritic 
values represent the initial values and are known with much better ac- 
curacy (often better than 10%) than the photospheric ones, we have 
assumed that the initial solar heavy metal abundances are given ap- 
proximately by the meteoritic (CI chondrites) values of Grevesse and 
Noels (1993a). 
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3.3.2 The CNO Abundances 



The careful and comprehensive analysis of atomic and molecular lines 
of C, N and O in the solar photosphere by Lambert (1978) has served 
as the standard source of their solar abundance. Improved atomic 
data, NLTE corrections and study of infrared lines have slightly changed 
the standard values of the photospheric CNO abundances to those 
given in Table I. These values are very much different from those found 
in primitive meteorites. In particular Carbon and Nitrogen are under- 
abundant in CI chondrites by about an order of magnitude compared 
with their photospheric abundances (e.g., Sturenberg and Howlwcgcr 
1990; Grevesse et al. 1990). We have assumed in our calculations that 
the initial CNO abundances are given by the photospheric abundances 
properly corrected for diffusion during the entire solar evolution. 

3.3.3 The Abundances of ^He and ^He 

Because the sun is a G2 star its photospheric Helium abundance can 
not be measured directly. The photospheric Helium abundance of 
much younger hotter B stars in the solar neighborhood can only be 
used for rough estimates of its initial solar abundance. The initial 
^He mass fraction in the solar nebula is known only approximately, 
0.24 < Y < 0.30. The predictions of the solar models are rather 
sensitive to the initial mass fraction of ^Hc. Consequently, the initial 
^He solar abundance has been treated in the standard solar models as 
an adjustable parameter. The present day ^He surface mass fraction 
can be inferred from helioseismology data which yield a surface ^He 
mass fraction of Yg = 0.242 it 0.003 (Hernandez and Christensen- 
Dalsgaard 1994). However, their formal error is highly misleading 
because of the great sensitivity of the result to the model of the solar 
atmosphere, the equation of state there and the atmospheric opacities. 
We estimate that at present the ^He mass fraction at the solar surface 
is not known from helioseismology better than Yg = 0.242 ± 0.030. 

Prom measurement of the [^He]/['^He] ratio in meteorites and in the 
solar wind and from the above rough estimate of the initial abundance 
of ^He, Geiss (1993) concluded that the presolar abundance of ^He is, 
pHe]/[H]=(1.5±0.3)xlO-5 (by numbers), i.e., log(pHe]/[H])+12=7.18± 
0.08. 
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3.3.4 The Deuterium Abundance 

Deuterium is easily destroyed already at relatively low temperatures. 
Consequently, all the primordial deuterium has been destroyed in the 
Hayashi phase and its initial abundance cannot be estimated from the 
solar photosphere. The Deuterium abundance in the local interstel- 
lar medium was recently measured with the Hubble Space Telescope 
(Linsky et al. 1993) to be, [D]/[H]=1.65l[J:?^ x 10"^ (by numbers), i.e., 
log([D]/[H])+12=7.22 ± 0.05. This value is consistent with the initial 
solar value, [D]/[H]=(2.6± 1.0) x 10-^ obtained by Geiss (1993) from 
the analysis of solar wind particles captured in foils exposed on the 
moon and from studies of primitive meteorites, which we used. 

3.3.5 The Abundances of ^Li, ^Be and ^^B 

The photospheric abundances of ^Li, ^Be and ^^B are smaller by a 
factor of nearly 150, 3 and 10, respectively, than their meteoritic 
abundances. The origin of such large differences is still not clear. 
However, the initial solar (meteoritic) abundances of Lithium, Beril- 
ium and Boron are very small and do not play any significant role in 
solar evolution. On the other hand, their depletion can provide a clue 
to the history of the convection zone and the solar evolution. 



4 The Physical Input 
4.1 Nuclear Reaction Rates 

The cross sections for most of the nuclear reactions that play impor- 
tant role in the sun fall steeply when the energy drops below the 
Coulomb barrier. At solar energies they become too small to be 
measured directly in laboratory experiments. Consequently, they are 
either calculated theoretically or extrapolated from laboratory mea- 
surements at much higher energies. Because the results of different 
experiments often differ considerably, and because different theoreti- 
cal models used for the extrapolation often yield different results, the 
nuclear reaction rates used in different standard solar models depend 
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on personal judgment. Below we discuss briefly some of the prob- 
lematics of obtaining reliable thermonuclear reaction rates for solar 
conditions and we explain our specific choices. These choices are also 
summarized in Table II and compared there with the choices of BP95. 
For all other reactions we have chosen to use the most detailed and 
complete compilation of thermonuclear reaction rates published by 
Caughlan and Fowler (1988). 
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4.1.1 The Problematics 



To eliminate from the extrapolation the strong energy dependence due 
to Coulomb barrier penetration, cross sections are normally parametrized 
as 

a.,{E) = M^e-^-'^ , (10) 

where 2iTr] = 31. 29 ZiZjy^fT/E is the Sommerfeld parameter, Zj and 
Zj being the charge numbers of the colliding nuclei, their reduced 
mass in atomic mass units and E their center-of-mass energy in keV. 
The factor Sij{E) is expected to vary rather slowly with energy. It 
is usually extracted from the measured cross section at laboratory 
energies and extrapolated to the lower solar energies, using either a 
polynomial fit or an energy dependence that follows from a specific 
model for the reaction. 

The uncertainties in the nuclear reaction rates at solar conditions 
are still large due to (1) uncertainties in the measured cross sections 
at laboratory energies (2) uncertainties in the extrapolations of the 
Sij{E) from laboratory energies down to solar energies, (3) uncertain- 
ties in dense plasma effects (screening, correlations and fluctuations) 
on reaction rates. Unfortunately, only for a few simple reactions ( 
p + p^D + e^ + t'e, p + e + p^D + fe and e~ + ''Be — > ''Li + v^) 
the cross sections can be calculated accurately from theory. For all 
other direct reactions neither the microscopic methods (for reviews see, 
e.g., Langanke 1991) such as the Resonating Group Method (RGM) 
and the Generator Coordinate Method (GCM), nor the potential mod- 
els such as the Optical Model (OM) and the Distorted Wave Born 
Approximation (DWBA), can predict accurate and reliable low en- 
ergy cross sections. For instance, the RGM and the GCM, which 
are currently considered to be the best theoretical methods for cal- 
culating direct nuclear reactions, predict (see, e.g. Descouvemont 
and Baye 1994, Johnson ct al. 1992) 517(0) f=i 25 - 30 eF • 6 for 
the reaction p -|- '''Be ^ ^B -|- 7. However a simple inspection of the 
results of these models reveals that they reproduce poorly the mag- 
nitude of the measured cross section, the position of the resonance, 
the width of the resonance, the height of the resonance and the ob- 
served shape of the cross section as function of energy. To avoid 
these discrepancies only the energy dependence of these models has 
been used by Johnson et al. (1992) to extrapolate the measured 
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cross sections to E = 0, yielding £'17(0) « 22.4 eV • b . This value 
has been used in BP92 and BP95. However, the procedure used 
by Johnson et al. (1992) is rather an ad hoc one and it is ques- 
tionable in view of the facts that (a) their model does not repro- 
duce accurately enough the measured energy dependence of the cross 
section at lab energies, and (b) they have extrapolated an "aver- 
age cross section" obtained by averaging cross sections which differ 
by many standard deviations and have different energy shapes: The 
cross sections measured by Kavanagh (1960) and by Parker (1968) 
differ by more than Scr from the cross sections measured later by 
Vaughn (1970) and by Filippone (1983) in the same energy ranges 
( [Kavanagh] /[Filippone] = 1.34 ± 0.11, [Parker] /[Vaughn] = 1.42 ± 0.13; 
see, e.g., Gai 1995). 

4.1.2 The Choice of Extrapolation 

Dar and Shaviv (1994) have pointed out that sub-Coulomb radiative 
captures and transfer reactions take place mainly when the colliding 
nuclei are far apart and not when their centers overlap. They argued 
that most of the energy dependence of the astrophysical Sij (0) factors 
is because Eq. 10 ignores this fact. The Coulomb barrier penetration 
factor exp{—2T:rj) / ^/E in Eq. 10 actually represents 

l^c(0)P-^^, (11) 

the absolute value squared of the Coulomb wave function at the ori- 
gin (normalized asymptotically to a plane wave of unit amplitude), 
while the physical cross section is proportional to the square of the 
Coulomb wave function at the effective distance R where the reaction 
takes place, i.e., to the Coulomb barrier penetration factor at a relative 
distance R. They concluded that the optical models which describe 
well the shapes of the bound state and relative motion wave functions 
outside the nuclear potential are the most reliable for extrapolating 
the laboratory cross sections to solar energies. Alternatively, they 
proposed that after extracting the energy dependence of the Coulomb 
barrier penetration factor to the effective distance R where the re- 
action takes place (and other trivial energy dependencies) a simple 
polynomial fit can be used to extrapolate the lab cross section to solar 
energies. 
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Let us show that the reactions take place mainly when the nuclei 
are far apart. Consider for instance radiative captures. The transition 
amplitude is proportional to an overlap integral I = J u* {r)ripc(f)dr 
where u and ipc arc the radial parts of the bound state wave function 
and of the Coulomb distorted wave function of the initial relative mo- 
tion, respectively. The bound state wave function outside the nuclear 
radius falls exponentially and has the asymptotic form, u{r) ~ e~^'^ 
where /3 = ^/2pEh/fi) with being the nuclear binding energy of 
the colliding nuclei in the final nucleus. The incident Coulomb wave 
function decreases inside the Coulomb barrier (in the WKB approxi- 
mation) like 

where Vc{r) is the effective Coulomb potential (Coulomb plus centrifu- 
gal) and Ro is the classical turning point. At keV energies, E is much 
smaller than Vc{r) near the nucleus and 

^^(r) ~ r+V4g-2»?[7r/2-ar-csm(E/Be)i/2-(E/-Ec)i/2(l-E/-Ec)i/2] _ rV4g7V^^ 

(13) 

where Ef. = ZiZje^ /r and 7 = JSZiZje^ji/h. Consequently, the con- 



tribution to the radial overlap integral comes mainly from the vicinity 
of r = i?, where 



As an example, for the reaction ^He +^ He Be -|- 7 one finds R ~ 

10 /m, while for the reaction p +^ Be -^^ B + 7 one finds i? « 70 fm 
in good agreement with the Optical Model and DWBA calculations. 
Eq. 14 is valid only if i?jv ^ R ^ Ro, where i?jv is the nuclear radius 
of the i-|-j bound state and Rq = ZiZj/E is the classical turning point. 

11 R < Rn then R ~ Rjsf. Similarly, ii R > Rq then R Rq. Thus, 
instead of being proportional to the Coulomb barrier penetration fac- 
tor at the origin, the radiative capture cross section is proportional 
to k^\u{R)Ripc{R)\'^ ■ One can divide out this energy dependence and 
use a polynomial fit to extrapolate the reduced cross section to solar 
energies (Dar and Shaviv 1994) or use instead the Optical Model or 
DWBA to extrapolate the total cross section to solar energies. The 
resulting astrophysical S'j.j factors are very similar. They are summa- 
rized here only for the main reactions. 
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4.1.3 SniO) 



Extrapolations of the cross sections measured by Vaughn (1970) and 
by Filipponc(1983), using either simple potential models or the above 
very general properties of sub-Coulomb cross section, gave 517(0) ^ 
17 eV ■ b (Barker and Spear 1986; Dar and Shaviv 1994; Kim et al. 
1994). This value seems to be supported also by other types of ex- 
periments: Analysis of recent measurements of the virtual reaction 
p+^Be through Coulomb dissociation of in the Coulomb 
field of Lead (Motobayashi et al. 1994) gave 517(0) 17 eF • 6. A 
similar value, 5*17(0) « 17.6 eV ■ b, was estimated from the virtual re- 
action pv -|-^ Be -^^ B measured through the proton transfer reaction 
3He+^Be^D+8B (Xu et al. 1994). Consequently, we a have adopted 
the value 5i7 = 17 eF • 6 in our standard solar model calculations. 
The value 5i7 = 22.4 eV ■ b was used in BP92 and BP95. 

4.1.4 ^34(0) 

The value 534(0) = 0.51±0.02 keV-b was obtained from measurements 
of the prompt 7-ray emitted in the reaction ^He He Be + 7, 
while measurements of the induced ^Be activity led to a weighted av- 
erage 534(0) = 0.58 ±0.02 keV ■ b, which are different by 3.5 standard 
deviations. The origin of this discrepancy is still not known (Hilge- 
meier et al. 1988). Normalization to known cross sections favor the 
lower value. 

Using the measured energy dependence of the cross section for 
^He +^ He — Be -I- 7 by Krawinkel et al. (1982) for extrapolating 
the cross sections measured by prompt gamma ray emission (Parker 
and Kavanagh 1963; Osborne et al. 1982; Krawinkel et al. 1982 
(multiplied by 1.4); Hilgemeier et al. 1988) to zero energy, we obtained 
534(0) = 0.45 keV ■ b. The value 534(0) = 0.533 keV ■ b was used in 
BP92 and 534(0) = 0.524 keV ■ b was used in BP95. 
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4.1.5 5'33(0) 



A similar analysis of the low energy data of Greife et al. (1994), 
Krauss et al. (1987) and Dawarakanath and Winkler (1971) on the 
reaction -^He He He + 2p gave 5*33(0) = 5.6 MeV-b. Essentially 
the same value, 533(0) = 5.57 MeV ■ b, was obtained by Krauss et al. 

(1987) and by Greife et al. (1944) by applying a polynomial fit to 
their data. The value ^33(0) = 5.0 MeV ■ b was used in BP92 and 
5-33(0) = 4.99 MeV ■ b was used in BP95. 

4.1.6 511(0) 

The cross section for the reaction p + p — > D + e+ + i^e is too small 
to be measured directly in the laboratory. The measured cross sec- 
tions for the weak isospin related reactions i/e+D— e^+n+n, i/e+D— 
'^e+P+n, and 7+D^n+p were used to obtain the relevant nuclear ma- 
trix element needed for calculating the cross section for p-|-p-^D-|-e"'" -|- 
fg. This procedure yielded a best value 5ii(0) 4.07 x 10~^^ keV ■ 
b which is consistent with the value used by Caughlan and Fowler 

(1988) . It is 4.2% larger than the value 5ii(0) ^ 3.896 x lO'^^ keV ■ b 
calculated recently by Kamionkowski &; Bahcall (1994) and used in 
BP95. 

Table II compares the astrophysical Sij{0) factors used in this work 
(and also in DS94) and in BP95. 

4.1.7 Screening Enhancement of Reaction Rates 

Screening of target nuclei by their electrons is known to enhance sig- 
nificantly laboratory nuclear cross sections at very low energies (e.g., 
Engstler et al. 1988) although a complete theoretical understanding 
of the effect is still lacking (see, e.g., Shoppa et al. 1993; Rolfs 1994 
and references therein). Screening corrections to the nuclear reaction 
rates in the SSM are usually represented by an enhancement factor 

where At/ is the gain in electronic potential energy when an incident 
ion of charge Zje penetrates the charged cloud around an ion of charge 
ZiC, T is the plasma temperature and Rd is the Debye length, Rd = 
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(fcr/47re2 2^271^)^2. This screening enhancement of thermonuclear 
reaction rates near the center of the sun is quite considerable, being 
5%, for the pp and pep reactions, 20% for the ^He^He, ^He^He and 
p^Be reactions, and 30%, 35% and 40% for the p capture by the C, N 
and O isotopes, respectively. 

The change in the Coulomb barrier (or equivalently the gain in 
energy) due to screening is estimated from the Debye-Hxickel (DH) 
approximation for the screened potential around a static ion in an 
electrically neutral plasma {TiZn^ = 0, being the number density 
of particles of charge Ze with Z = —1 fov electrons): 

= ^^-r/Ro ^ _ ^ for ^ « i?^, (16) 

r r Rd 

which is an approximate solution to Poisson's equation 

V^$i = 47reSZn^ = 47reSZn^e-^^*'/'=^ fii Rjj^i, (17) 

near a static ion of charge Zj. This equation and solution, however, are 
valid only far from the nucleus where eZ^i <C kT, or else the expan- 
sion is not valid. In the core of the sun where kT (~ 1 keV) is much 
smaller than the Coulomb barriers 1 MeV) between the reacting 
nuclei, most of the contribution to the nuclear reaction rates comes 
from collisions with cm. energies E S> kT. At the classical turning 
point eZj^i = E. Consequently, inside the barrier, eZj^i ^ kT, and 
the naive use of the DH approximate potential for calculating the bar- 
rier penetration factor (e.g. Clayton 1968) is unjustified for projectiles 
with wavelengths much shorter than the classical distance of closest 
approach. The DH solution also assumes that there are many ions 
and many electrons within a Debye sphere. Also these conditions are 
not quite satisfied in the core of the Sun, where the inter-ion distance 
is about the same as the Debye length. It was found that when the 
conditions for the validity of the DH approximation arc violated in 
laboratory plasmas, calculations based on the DH approximation fail 
dramatically (Goldsmith et al. 1984) in reproducing various obser- 
vations, unless the wavelength of the particle under consideration is 
much larger than the Debye length. The later is the case for the elec- 
tron capture reaction ^Be(e~, fe)^Li, where it was found by Johnson 
et al. (1992) that the DH screening potential changes |V'e(0)P by only 
a few percent. However, for the thermonuclear reactions in the sun the 
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wave lengths of the reacting ions are very small in comparison with 
their classical distances of closest approach and the Debye length. 

Moreover, an ion approaching from infinity a target nucleus gains 
from the Debye cloud a potential energy MJ e^ZiZj/Ro when it 
penetrates it. But, there is no gain in potential energy if either its 
initial position is already inside the Debye cloud (where the cloud 
potential is constant V eZi/Rr)) or if the ion leaves and enters 
similar potential wells. 

Near the center of the Sun, where pc ~ 156 g cm~^, Tc ~ 1.57 x 
lO^K, XciH) « 34%, and Xc{He) « 64%, the average inter-ionic 
spacing is n~^^^ ~ 2.8 x 10~^ cm, similar to the Debye length, Rd ~ 
2.3 X 10~^ cm. On the other hand the gain in potential energy is in- 
creased if both ions retain their Debye clouds during their approach. 
However, the most effective ion energies are 5-20 kT and their veloci- 
ties are not much smaller than the average thermal electron velocities. 
There may not be sufficient time for the plasma to rearrange itself 
around the fast moving ions and screen them effectively (the typical 
rearrangement time, 27r/LUp, where cOp = (47rng/me)^/^, is longer than 
the time it takes a fast ion to cross the Debye length) . All these effects 
may modify substantially the screening enhancement of the fusion re- 
action rates near the center of the Sun as calculated from the static 
DH screening potential. Actually, the situation is too complicated for 
reliable analytic estimates. Reliable estimates, however, may be based 
on numerical N body simulation of the classical trajectories of elec- 
trons and ions in a dense plasma which can be used to evaluate the 
various effects of correlations and fluctuations on the thermonuclear 
reaction rates in dense plasmas (Shaviv and Shaviv 1995). To test the 

sensitivity of the standard solar models to the screening corrections 
we have carried out the calculations with and without the standard 
screening enhancement factors of all the thermonuclear reaction rates. 
We found (Dar and Shaviv 1994) that removing/including the screen- 
ing enhancement factors for all nuclear reaction rates had only a small 
effect on the calculated solar neutrino fluxes, due to accidental cancel- 
lations. Screening enhancement may, however, play important role in 
other stars. Moreover, changes in the screening factors which change 
their ratios for the different reactions may cause significant changes in 
the calculated solar neutrino fluxes. 
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4.1.8 Radiative Opacities 



The radiative opacities depend on the local chemical composition, den- 
sity and temperature in the sun. We have used radiative opacity tables 
computed by the OPAL group at Lawrence Livermore National Lab- 
oratory that were kindly provided to us by F.J. Rogers. They are 
updated version of the OPAL tables of Rogers and Iglesias (1992) 
for the most recent determination by Grevesse and Noels (1993a) of 
the heavy element composition of the sun from the meteoritic and 
photospheric data. As for the low temperature opacities, we used the 
Alexander, Rympa & Johnson (1983) low temperature opacities which 
include the effect of molecules. Note that the original table contains 
a numerical error. This error was corrected by interpolation from ad- 
jacent values. Since the opacities are differentiated for the numerical 
scheme a spline interpolation is used to guarantee regular derivatives 
and opacities. Note that the nuclear reactions and the element diffu- 
sion change the heavy element abundance throughout the sun as well 
as its internal breakdown to the various elements. The breakdown 
of the OPAL opacities to the individual contributions of the various 
elements as function of local density, temperature and composition is 
not available yet for general use. However, the calculations of Kovetz 
and Shaviv (1994) show that the diffusion of the heavy elements is 
approximately the same for all the heavy metals including the partial 
ionization of each element (see e.g.. Table I of Kovetz &: Shaviv 1994). 
Thus, the error introduced in the radiative opacity by the neglect of 
the change in the heavy element breakdown is expected to be negligi- 
ble compared to other inaccuracies in the opacities. Besides diffusion 
the only important process that changes the relative abundances of the 
heavy elements is the conversion of ^^C to ^^N in the core of the sun. 
However, since the radiative opacities of C and N are not very different 
when they are completely ionized, we do not expect this conversion to 
affect significantly the total radiative opacity. Consequently, we cal- 
culated the opacities for the values of p, T, X and Z in each shell by 
interpolating the OPAL opacities, assuming the same breakdown of Z 
to the heavy elements. We plan to extend the present calculations to 
include the effects of the variations of the heavy element breakdown 
on the radiative opacity. 
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4.1.9 Equation of State and Partial Ionization 



Under the conditions prevailing in the core of the sun, the electrons 
arc partially degenerate. Their contributions to the various thermody- 
namic quantities (pressure, entropy, adiabatic exponents) is calculated 
in the Kovetz Shaviv code by exact evaluation of the Fermi-Dirac 
integrals. The Coulomb interactions in the solar plasma have been 
included in the non ideal part of the free energy. The value of the 
Coulomb parameter in the core of the sun is about 0.1-0.2. Kovetz & 
Shaviv (1994) and Dar & Shaviv (1994) used a fit by Iben et al. (1992) 
to the results of Slattery, Doolen Sz De Witt (1980) for the internal 
energy of the gas which were obtained for the Coulomb parameter 
r < 1. These results were calculated for a One Component Plasma 
(OCP) and ignore the contribution of the electrons. Recent results 
by Shaviv and Shaviv (1996) indicate that the assumption of OCP is 
not so good for such low values of F and hence, we calculated a model 
with the simple Debye-Hiickel model. The local ionization state of 
each element is calculated by solving the local Saha equations. Par- 
tial ionization effects are included in the equation of state and in the 
diffusion coefficients (C,N, O or Mg are once to twice ionized at the 
surface and completely ionized in the core. Fe is doubly ionized at 
the surface but only 16-20 times ionized in the core). Partial ioniza- 
tion has an important effect on the equation of state in the outer part 
of the sun even if the relevant region is convective since the partial 
ionization affects the entropy of the gas. 



5 The Solar Evolution Code 



Modeling the solar evolution requires integration of the standard stel- 
lar evolution equations over very long times. Truncations and inaccu- 
racies in the time integration may add to large errors if not properly 
handled. In particular, the numerical integration must ensure energy 
conservation. Namely, the integral of the total solar luminosity in 
photons and neutrinos over a long time must be equal to the change 
in the total solar binding energy (nuclear, gravitational, thermal and 
Coulomb) during that time. A carefully constructed advanced stellar 
evolution code was described by Kovetz and Shaviv (1994). We have 
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used this code to calculate the evolution of a 1 Mq protostar, fully con- 
vective with a uniform composition, initial solar element abundances 
and a luminosity around 12Lq. The code follows the Hayashi track, 
the settling to the main sequence (ZAMS) and the main sequence 
evolution. It includes hydrodynamic, gravitational and nuclear en- 
ergy release, radiative energy transport, convection, and diffusion of 
all elements more abundant than a prescribed abundance. Neither 
nuclear equilibrium, nor complete ionization are imposed. The effects 
of partial ionization are included in the equation of state, the diffusion 
coefficients and the radiative opacities of OPAL and in the calculation 
of the solar atmosphere. For the convenience of the reader we high- 
light here some important details of the stellar evolution code. For 
more details the reader is referred to its description in the original 
paper of Kovetz and Shaviv (1994). 
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5.1 Nuclear Evolution 



The nuclear evolution equations for the local number fractions yi = 
Xi/Ai have the general form 



where S{y) is a matrix with positive (creation) and negative (destruc- 
tion) elements. Most of these elements are linear in the y's, but some 
(corresponding to beta decays) are constants. Thus, the right hand 
side of Eq. 18 is quadratic in the y's. The basic nuclear network con- 
tains all the isotopes and nuclear reactions which appear in Caughlan 
and Fowler (1988). The nuclear reactions rates have been updated 
with the reaction rates of Wiescher et al. (1989), Kubono et al. (1989) 
and Gorrcs ct al. (1989). The other changes (updating the normaliza- 
tion of the '^Hc'^Hc, '^Hc'^Hc, p^Be reactions) were described in section 
4. The equation is solved by conversion into a finite difference equa- 
tion. Stability and the requirement that the solution approaches the 
solution for nuclear equilibrium as time tends to infinity require that 
the time differencing be implicit. Thus, Eq. 18 is replaced by the fully 
implicit finite difference system. 



which is correct to 0{6t), numerically stable and delivers the correct 
values yi{t + St) for those isotopes that change very rapidly and are 
thus in nuclear equilibrium. No special provision is required for han- 
dling isotopes that are in or tend towards nuclear equilibrium. The 
non linear equations 19 are solved by iterations in which large ma- 
trices are inverted (which consume a lot of CPU time). About three 
iterations were needed to obtain an accuracy of 1 : 10^^*^ needed in 
order to follow the concentration of species like ^Li and ®B. Our nu- 
clear network contains about 80 isotopes, of which only the following 
18 isotopes were found to play important role in the nuclear evolution: 



H,2 D,3 He,4 He,7 Li,^ Be,^ B,'^ C,'^ C,''' N,^^ N,!^ N,!^ 0,1^ 0,1^ 0,1^ F,'^ F,20 Ne, 



^Na. The nuclear network searches automatically for the most impor- 
tant contributing reactions at each phase. In the particular case of the 
sun the network was reduced to these 18 isotopes. Yet, reactions like 
^He^D were included by the program. (The other elements, in par- 
ticular Ne, Na, Mg, Si, S and Fe, which are present in the sun in 



(18) 



yt+8t = yt + StS{yt+st)yt+5t, 



(19) 
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significant quantities that affect the "heavy element" abundance Z 
but do not participate in the nuclear reactions, were included in the 
calculations of the opacities, the equation of state, the thermodynamic 
quantities and the diffusion of elements). The nuclear energy genera- 
tion rate Qn^ which is the time derivative of the total nuclear binding 
energy is calculated from: 

Qn = -EbS(yt+5t)yt+5t, (20) 

where is the vector of nuclear binding energies. This form ensures 
that the integral over time of the nuclear energy released in the sun 
is equal to the change in the solar nuclear binding energy during that 
time. 

5.2 Space and Time Steps 

The sun was divided into 550 mass shells, distributed in such away 
to yield sufficiently fine distribution in the core, near the bottom of 
the outer convective zone (where Am was 1 x 1Q~^Mq) and at the 
surface. The time steps were constrained in such a way that the 

changes in logT, logp and logXi during a time step were smaller than 
0.02, 0.05 and 0.02, respectively. No time steps larger than 3 x 10^ y 
were allowed. The implicit difference equations, replacing the laws of 
energy and momentum conservation were iterated until logT and logp 
were determined within 1 x 10~^. 

5.3 Convection 

A region in which dp /dp < 0, or, in the case of uniform composition 
and ionization, a region where the entropy decreased outwards, was 
considered to be convective. In a convective region a complete and 
instantaneous mixing was carried out and a convective flux was added 
to the radiation flux. The convective flux was calculated according 
to the mixing length recipe, using the numerical constant of Mihalas 
(1973). 
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5.4 The Surface Boundary Condition 

The outer half of the last mass shell was treated as a static atmosphere. 
The Stefan-Boltzmann luminosity from the last shell was equated to 
the total heat flow into the shell and the effective temperature was 
adjusted (by iterations) until the calculated entropy (as given by the 
model atmosphere) at the base of the atmosphere was equal to the 
entropy (as given by the evolution program) at the center of the shell. 
The Burlish-Stoer extrapolation method (Press et al. 1989) was used 
for the inward integration. A fine grid of model stellar atmospheres is 
constructed and the conditions at the middle of the last mesh point 
are found by interpolation inside the this grid. The model atmosphere 
start with the photosphere r = 2/3 and no T — r relation is required. 
The equations of the gray atmosphere are integrated down to the mass 
of the prescribed (half) last mass shell. 

5.5 Diffusion 

Diffusion, caused by density, temperature, pressure, chemical compo- 
sition and gravitational potential gradients may play an important 
role in the sun. The effects of helium diffusion on solar models were 
first studied by Noerdeling (1977). Improved calculations were done 
by Cox et al. (1989), ProfHtt and Michaud (1991), and by Bahcall 
and Pinsonncault (1992) who incorporated the approximate analyti- 
cal description of hydrogen and helium diffusion by Bahcall & Loeb 
(1990) in their standard solar model calculations. Proffit (1994) and 
Kovetz & Shaviv (1994) have studied the effects of diffusion on the 
sun by solving numerically the diffusion equations for all the elements 
with mass fraction greater than a prescribed value (usually 10~^). To 
simplify the numerical handling of the ionization, each isotope was as- 
sumed to be in a single ionization stage at each point, the one which is 
the most abundant. Consequently, the ionization stage changed from 
one shell to the other but in each shell only one state was allowed. 
The Kovetz-Shaviv code calculates the diffusion of all the individual 
elements from the premain sequence phase to the present age. The 
binary and thermal diffusion coefficients depend on the squared ionic 
charges. The ionization state of each element in every shell is cal- 
culated by solving the Saha equations for all the elements in each 
shell. All elements with mass fractions less than 10~^ are treated as 
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trace elements. By trace element we mean that collisions with other 
trace elements are neglected. As a consequence of diffusion the surface 

and internal element abundances change for each element in a slightly 
different way. Diffusion depletes the surface abundances of ^He and 
the heavy elements. For C, N and O only their present photosphcric 
abundances are known (see chapter 3.3.2) Consequently, in our solar 
model calculations their initial solar abundances were adjusted to re- 
produce their observed photospheric abundances. Thus, their initial 
abundances should be considered as a prediction of the standard solar 
model. For all other elements the initial meteoritic abundances are 
used to predict their final surface abundances which can be compared 
with their observed photospheric abundances. The unknown initial 
abundance of '^He is treated as a free parameter. It is adjusted to best 
reproduce the presently observed sun. Its predicted surface abundance 
today can be compared with the value derived from helioseismology 
(Hernandez and Christensen - Dalsgaard 1994). 



6 Results and Comparison With Ob- 
servations 

Our predictions of the solar neutrino fluxes for three standard so- 
lar models are summarized in Table Ilia and are compared with the 
results from the four solar neutrino experiments. The models were cal- 
culated with the stellar evolution code of Kovetz and Shaviv (1994) as 
described in section 5, with the improved astrophysical and physical 
input data which were described in sections 3 and 4. The three mod- 
els differ only in their treatment of element diffusion and equation of 
state: The model labeled DSND does not include element diffusion. 
The models labeled DS94 and DS95 include element diffusion while 
the model DS95 includes an improved equation of state for F < 1. 
Model DS94 assumes that the initial heavy metal abundances in the 
sun were equal to their meteoritic values, and those of C, N and O 
were equal to their observed photospheric abundances, as summarized 
in Table I. In model DS95 the initial abundances of C, N and O were 
adjusted to yield their present day photospheric abundances while the 
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heavy metal abundances were assumed to be those found in primi- 
tive meteorites, as summarized in Table I. Table Illb presents some 
physical characteristics of the three models. The indexes have the 
following meanings: c means the core, s the surface and means the 
initial values. Rconv is the radius of the base of the convective enve- 
lope and Tfjc is the temperature at the base of the convective envelope. 
As can be seen from Table Ilia, all three models yield a solar neu- 
trino flux consistent with that measured by Kamiokande. However 
all three models predict capture rates in the Chlorine and Gallium 
experiments which are significantly larger than those measured by 
Homestake, GALLEX and SAGE (We have used the neutrino cross 
sections from Table 8.2 of Bahcall (1989) to convert solar neutrino 
fluxes to capture rates). Tables IV present a comparison between four 

solar models. The model labeled BP92 is the best model of BP92. It 
includes diffusion of protons and ^He but not of other elements. The 
model labeled BP95 is the best model of BP95 which includes also 
diffusion of the heavy metals but assumes that all the heavy elements 
diffuse like fully ionized iron. The predictions of the DS models differ 
significantly from the BP models because they differ in input physics, 
approximations and numerical methods. Most of the differences are 
due to the use of different reaction rates as summarized in Table II, 
the different treatments of diffusion and the equation of state. This 
is demonstrated in Table V where we present a comparison between 
the best model of BP95 without diffusion labeled BP95ND, and a 
solar model, labeled DS(BPND) calculated with the Kovetz- Shaviv 
stellar evolution code with the same physical and astrophysical input 
parameters and the same nuclear reaction rates used in BP95ND and 
without inclusion of element diffusion. As can be seen from Tables Va 
and Vb the two calculations yield similar results. Even the fluxes of 
^B and CNO solar neutrinos which, under the imposed solar boundary 
conditions, are very sensitive to the central solar temperature differ by 
less than 4%. These remaining differences are mainly due to the use of 
different equations of state, numerical methods, fine zoning and time 
steps in the two codes and due to the inclusion of premain sequence 
evolution in our code. To emphasize the important role that might be 

played by diffusion. Tables Va and Vb also include the current best 
solar models of Dar and Shaviv (1995) and of Bahcall and Pinson- 
neault (1995), which include diffusion of all elements. As can be seen 
from Tables Va, Vb, Bahcall and Pinsonneault (1995) found rather 
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large increases in their predicted Be, B, ^N, and F solar neu- 
trino fluxes; 14%, 36%, 52%, 58%, and 61%, respectively, compared 
with their model (BP95ND) with no diffusion. These induce 36%, 
33%, 9% increases in the predicted rates in Kamiokande, Homestake, 
and in GALLEX and SAGE, respectively. However, we predict more 
moderate increases due to diffusion, 4%, 10%, 23%, 24% and 25%, 
respectively, in the above fluxes, which correspond to 10%, 10% and 
2% increases in the predicted rates in Kamiokande, Homestake, and 
in GALLEX and SAGE, respectively. The differences in the effects of 
diffusion in DS94 and BP95 are mainly due to two reasons: (a) In the 
calculations of Bahcall and Pinsonneault (BP95) all heavy elements 
were assumed to diffuse at the same rate as fully ionized iron while 
the Dar-Shaviv calculations (DS95) followed the diffusion of all the 
elements separately and used diffusion coefficients for the actual ion- 
ization state of each element, (b) Bahcall and Pinsonneault assumed 
that the meteoritic abundances represent the solar surface abundances 
today and not their initial values. They adjusted their initial values 
to reproduce surface abundances today equal to the meteoritic values. 
Thus they have actually used an initial ratio ZjX = 0.0285 (see Table 
IV of BP95) while the observed ratio in meteorites is ZjX = 0.0245 
(Grevesse and Noels 1993). Dar and Shaviv (1994; 1995) used the me- 
teoritic values for the initial metallic abundances and predicted present 
day depleted surface abundances. Unfortunately, the uncertainties in 
the measured photospheric CNO abundances (typically 30%-40%) arc 
much larger than their predicted depletion (typically 10%) and do not 
allow a reliable test of this prediction. 

The photospheric abundances of ^Li, ^Be and B are much smaller 
than their abundances in primitive meteorites. For instance, the pho- 
tospheric abundance of '''Li is smaller by nearly a factor of 150 (An- 
ders and Grevesse 1989). These differences are not explained by the 
standard solar models (see e.g., Morrison 1992), in spite of the fact 
that significant Lithium burning takes place during the Hayashi phase. 
However, observations of ^L in young solar mass stars (e.g., Soderblom 
et al. 1993) and solar evolution calculations (e.g. Kovetz and Shaviv 
1994) suggest that premain sequence destruction of ^Li together with 
its main sequence destruction is less than a factor 5 and cannot explain 
a surface depletion of ^Li by nearly a factor 150. Its destruction (via 
^Li-hp^^He-h^He) requires higher temperatures than ~ 2.2 x W^K 
predicted by the standard solar models for the bottom of the con- 
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vective zone. Differential rotation below the base of the convective 
zone may cause mixing and bring ^Li much deeper where it can burn, 
(see, e.g., Vauclair and Charbonnel). Such mixing may , however, in- 
hibit the inward diffusion of the heavy elements. Thus, it is not clear 
whether solar models with diffusion (such as DS95) provide a more 
realistic description of the sun than solar models without diffusion 
(such as DSND). Helioseismology data, however, is better explained 
by SSM's with diffusion. See, e.g., Cristensen-Dalsgaard et al. 1993. 
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7 Where Are The Be Neutrinos ? 

As we have seen in chapter 6, standard solar models predict cap- 
ture rates in the Chlorine and Gallium experiments which are sig- 
nificantly larger than those measured by Homcstakc, GALLEX and 
SAGE. Combined with the results from Kamiokande, they seem to 
suggest that the ^Be solar neutrino flux is strongly suppressed (and 
perhaps also the fluxes of and CNO solar neutrinos). Can standard 
physics explain the suppression of the '^Bc solar neutrino signal in ^^Cl 
and ^^Ga ? We think that such a possibility has not been ruled out 
yet: Electron capture by ^Be into the ground state of ^Li produces 
862 kcV neutrinos. The threshold energy for neutrino absorption by 
^^Cl is 814 keV. Thus, absorption of ^Be neutrinos by ^^Cl produces 
48 KeV electrons. The average energy of the pp solar neutrinos is 
265 KeV. The threshold energy for neutrino absorption in ^^Ga is 233 
KeV. Consequently, the produced electron has a typical energy of 33 
keV. The de Broglie wave lengths of such electrons are larger than 
the Bohr radii of the atomic K shells in CI and Ga and their ener- 
gies are similar to the kinetic energies of electrons in the K shells. 
Consequently, screening of the nuclear charge by atomic electrons and 
final state interactions (exchange effects, radiative corrections, nuclear 
recoil against the electronic cloud, etc.) may reduce the absorption 
cross sections of pp neutrinos in ^-"^Ga (perhaps making room for the 
expected contribution of '^Be in Gallium ?) and of ^Bc neutrinos in 
^^Cl (perhaps making the solar neutrino observations of Kamiokande 
and the Homestake experiment compatible Dar, Shaviv and Shaviv 
1996). It is interesting to note that although final state interactions in 
Tritium beta decay have been studied extensively, they do not explain 
yet why its end-point spectrum (E^ 18.6 kEV) yields consistently, 
in all recent measurements, a negative value for the squared mass of 
the electron neutrino. Final state interactions in ^^Cl and ^^Ga are 
expected to be much larger because of their much larger values of Z. 
Note also that this explanation of the solar neutrino problem implies 
that experiments such as BOREXINO and HELLAZ will observe the 
full ^Be solar neutrino flux while the MSW solution predicts that it 
will be strongly suppressed. 

Even if the ^Be solar neutrino flux is strongly suppressed, it does 
not eliminate yet standard physics solutions to the solar neutrino 
problem. For instance, collective plasma physics effects, such as very 
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strong magnetic or electric fields near the center of the sun, may po- 
larize the plasma electrons, and affect the branching ratios of elec- 
tron capture by ^Be (spin 3/2~) into the ground state (spin 3/2~, 
E^,^ = 0.863 MeV, BR=90% and the excited state (spin 1/2", E^,^ = 
0.381 MeV, BR=10%) of ^Li. Since solar neutrinos with Ey^ = 
0.381 MeV are below the threshold (0.81 MeV) for capture in ^^Cl 
and have a capture cross section in ^^Ga that is smaller by about a 
factor of 6 relative to solar neutrinos with Ey^ = 0.863 MeV , there- 
fore a large suppression in the branching ratio to the ground state can 
produce large suppressions of the ^Be solar neutrino signals in ^^Cl 
and in ^^Ga. 



8 The MSW Solution 



Standard solar models, like the one presented in this work, perhaps 
can explain the results reported by Kamiokande. However, standard 
physics cannot explain an ^^Ar production rate in ^^Cl smaller than 
that expected from the solar ^B neutrino flux measured by Kamiokande 
If the experimental results of Kamiokande and Homestake are inter- 
preted as an evidence for such a situation (e.g., Bahcall 1994; 1995), 
they do imply new physics beyond the standard particle physics model 
(Bahcall and Bethe 1991). In that case an elegant solution to the 
solar neutrino anomaly is resonant neutrino flavor conversion in the 
sun, first proposed by Mikheyev and Smirnov (1986) (see also Wolfen- 
stein 1978; 1979). Many authors have carried out extensive calcula- 
tions to determine the neutrino mixing parameters which can bridge 
between the predictions of the standard solar models and the solar 
neutrino observations. They found that a neutrino mass difference 
Am^ ~ 0.7 X 10~^ eV"^ and a neutrino mixing of siv?29 k, 0.5 x 10~^ 
can solve the solar neutrino problem (see, e.g., Krastev and Petcov 
1993; Castellani et al. 1993; Hata et al. 1994; Kwong and Rosen 
1994; Berezinsky et al. 1994). In fact, these parameters can be found 
from simple analytical considerations: As we have seen, the results 

of GALLEX and SAGE can be predicted directly from the pp flux, 
which is fixed by the solar luminosity, and the ^B flux measured by 
Kamiokande, if GALLEX and SAGE are blind to all other neutrinos. 



33 



This blindness may be due to flavor conversion of these neutrinos. The 
resonance condition for flavor conversion in the sun is 

ne = rf/^ = — ^ — ^ , (21) 

for small mixing angles. Solar pp neutrinos, whose maximum energy is 
0.42 MeV do not encounter resonance density if the resonance density 
for Ey = 0.42 is larger than the central electron density in the sun. 
The 0.861 MeV ^Be solar neutrinos, and all the other solar neutrinos 
which are more energetic, encounter a resonance density in the sun 
and suffer a resonance flavor conversion if the resonance condition is 
already satisfied for 0.861 MeV neutrinos. The electron density at the 
center of the sun IS Tip 6 X 10^^ cm *^ and the last two conditions 
together with Eq. 21 yield 

0.5 X 10^ eV'^ < Am^ < 1 x 10^ eV'^. (22) 

The probability of flavor conservation is given approximately by (e.g. 
Parke 1986, Dar et al. 1987) P{ue ^ u.^) e""/-^-, where 

e = , (23) 

where H is the solar scale height at the resonance. Thus, we see 
that the suppression of the solar Vg flux can still be adjusted by the 
choice of the mixing angle 9. A suppression factor of P at 
the most effective solar neutrino energies (~ 10 MeV) in the chlorine 
detector yields a suppression factor of e"^'^ at E^ ^ 1 MeV. Such a 
suppression and Am^ Ri 0.75 x 10~^ eV"^ yield a mixing parameter 
sin'^2e Ri 0.5 X 10-^. 

Are neutrino oscillations responsible for the solar neutrino anomaly? 
The answer may be provided by two experiments planned to begin col- 
lecting data in 1966: The Sudbury Neutrino Observatory (SNO) will 

detect two neutrino interactions in 1000 tons of heavy water: 

Ue + D^p + p + e- , (24) 

and 

Vx+D^p + n + Ux ■ (25) 

If more solar neutrinos are detected via the second reaction (which 
is blind to the neutrino flavor) than via the first reaction (where only 
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VeS can contribute), it would provide evidence that some solar i/^s 
have changed their flavor before reaching the detector. The total solar 
neutrino flux will also be measured at SNO by elastic scattering on 

electrons in the light water shield. The expected counting rate in 
the 1000 tones of heavy water of SNO is about 10 per day!. The 

Superkamiokande experiment will utilize 50,000 tones of high purity 
light water and will enhance the counting rate of solar neutrinos in 
Kaniiokandc by about a factor 30. It should have enough statistics 
to detect the change in the energy spectrum of the high energy 
solar neutrinos predicted by the MSW (Mikheyev and Smirnov 1986, 
Wolfenstein 1978) solution to the solar neutrino problem. 



9 Summary and Conclusions 



The results of the four pioneering solar neutrino experiments, confirm 
that the sun derives its energy by fusion of hydrogen into helium in 
its core. This is a great triumph both for experimental and theo- 
retical physics. The reliability of the results from these difficult and 
ingenious experiments is supported by their consistency, and recently, 
most convincingly, by the GALLEX Chromium source experiment . 
The capture rates of solar neutrinos measured by GALLEX and SAGE 

are above the minimal rates expected from the solar luminosity and 
the conservation of lepton flavor. The solar neutrino flux pre- 
dicted by our improved standard solar model, which was described 
in this paper, is consistent, within the theoretical and experimental 
uncertainties, with the solar neutrino observations at Homestake and 
Kamiokande. However, the experimental results from the four solar 
neutrino experiments seem to suggest that the flux of ^Be solar neutri- 
nos at Earth is much smaller than that predicted by the standard solar 
models, including our improved SSM which includes premain sequence 
evolution, element diffusion, partial ionization effects, all possible nu- 
clear reactions between the main elements, uses updated values for the 
initial solar element abundances, the solar age, the solar luminosity, 
the nuclear reaction rates and the radiative opacities and does not 
impose either nuclear equilibrium or complete ionization. Solutions 
to the solar neutrino problem which do not invoke neutrino properties 
beyond the standard electroweak model are not ruled out yet: 
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The solar neutrino problem may be a terrestrial problem. The 
neutrino capture cross sections near threshold in the radiochemical 
experiments may be different from the calculated cross sections (Dar, 
Shaviv and Shaviv 1996). The inferred solar neutrino fluxes from the 
GALLEX and HOMESTAKE experiments may be different from the 
true solar neutrino fluxes. They do not establish beyond doubt that 
there is a real ^Be solar neutrino deficit. Perhaps, future experiments 
such as BOREXINO and HELLAZ will be able to establish that. 

The solar neutrino problem may be an astrophysical problem. The 
deviations of the experimental results from those predicted by the 
standard solar models may reflect the approximate nature of the stan- 
dard solar models (which neglect, or treat only approximately, many 
effects and do not explain yet solar activity nor the surface depletion 
of Lithium, Berilium and Boron relative to their meteoritic values, 
which may or may not be relevant to the solar neutrino problem). 
Improvements of the standard solar model should continue. In par- 
ticular, dense plasma effects on nuclear reaction rates and radiative 
opacities, which are not well understood, may strongly affect the SSM 
predictions and should be further studied, both theoretically and ex- 
perimentally. Relevant information may be obtained from studies of 
thermonuclear plasmas in inertial confinement experiments. Useful 
information may also be obtained from improved data on screening 
effects in low energy nuclear cross sections of ions, atomic beams and 
molecular beams incident on a variety of gas, solid and plasma targets. 

Better knowledge of low energy nuclear cross sections is badly 
needed. Measurement of crucial low energy nuclear cross sections by 
new methods, such as measurements of the cross sections for the ra- 
diative captures p Be -^^ B + 7 and ^He +^ He Be + 7 by pho- 
todissociation of ^B and ^Be in the coulomb field of heavy nuclei are 
highly desirable. 

The ^^Ar production rate in ^^Cl may indeed be smaller than that 
expected from the flux of standard solar neutrinos as measured by 
electron scattering in the Kamiokande experiment. In that case neu- 
trino oscillations, and in particular the MSW effect, may be the correct 
solution to the solar neutrino problem. Only future experiments, such 
as SNO, Superkamiokande, BOREXINO and HELLAZ, will be able 
to supply a definite proof that Nature has made use of this beautiful 
effect. 
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Table I: Summary of Information on Abundances of Various Ele- 
ments relative to Hydrogen ( A=log([A]/[H])+12) in Primitive Mete- 
orites, in the Solar Photosphere, in the Solar Wind and in the Local 
Interstellar Medium, Used In The DS Standard Solar Models. 
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Element 


Abundance 


Source 


Reference 


D 


7.22 ± 0.05 


Meteorites, SolarWind 


Linsky 1993, Geiss 1993 


TT 

•^He 


7.18 ± 0.08 


Meteorites, Solar Wind 


Geiss 1993 


7-1- . 


1.54 ± O.OX 


Meteorites 


Anders and Grevesse 1989 


^Be 


1.13 ± O.OX 


Meteorites 


Anders and Grevesse 1989 




8.55 ± 0.05 


Photosphere 


Grevesse and Noels 1993 




6.60 ± 0.05 


Photosphere 


Grevesse and Noels 1993 


1 /I H. -T- 


7.97 ±0.07 


Photosphere 


Grevesse and Noels 1993 


IbQ 


8.78 ±0.07 


Photosphere 


Grevesse and Noels 1993 


9n-i*.T 

"'"Ne 


8.08 ± 0.06 


Photosphere 


Grevesse and Noels 1993 




6.33 ± 0.03 


Meteorites and Photosphere 


Grevesse and Noels 1993 


24]V[g 


7.58 ± 0.05 








6.47 ±0.07 






9Q j-^ . 


7.66 ± 0.05 






OlTD 


5.45 ± 0.04 








7.21 ± 0.06 








5.5 ±0.3 








6.52 ± 0.1 






40 


6.36 ± 0.02 






40k 


4.85 








3.08 






48Ti 


5.02 ± 0.06 






50v 


3.99 






52 


5.67 ± 0.03 






55 Mn 


5.39 ± 0.03 






56pg 


7.50 ± 0.04 








4.15 






58Ni 


6.25 ± 0.04 








4.33 






z/x 


0.0245 







Table II: Comparison Between the Astrophysical S Factors for 
the pp-chain Reactions used in BP95 and in DS94 and DS95. The 
values of S are given in keV ■ b Units. 
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Reaction S^P(O) S°^(0) 

^H(p,e+z/c)D 3.896 x lO'^^ 4.07 x lO"' 

^H(pe-z/e)D Bahcall 89 CF88 

3He(3He, 2p)He^ 4.99 x 10^ 5.6 x 10^ 

^Hej^He, 7)^Be 0.524 0.45 

^Bc(c-,z^e)Li'^ Bahcall 89 CF88 

^Be(p,7)8B 0.0224 0.017 
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Table Ilia: Comparison Between Solar Neutrino Fluxes Pre- 
dicted by the current best Standard Solar Models of Dar and Shaviv, 
with and without element diffusion, and the Solar Neutrino Observa- 
tions. 



V Flux 






DSND 


DS94 








6.10 


6.06 


(t)y(pep) [lO^cm" 






1.43 


1.42 








4.03 


4.00 








2.54 


2.60 




-^s 




3.21 


3.30 




-2s- 




3.13 


3.19 




-2s- 




3.77 


3.84 


[SNU] 






4.2 ± 1.2 


4.2 ± 


E(0a)Ga [SNU] 






116 ±6 


116 ± 


S(0(7)Ga [SNU] 






116 ±6 


116 ± 



DS95 Obs Exp 

6.10 

1.43 

3.71 -^(f^f^CBe) ALL 

2.49 2.9 ± 0.4 Kam. 

3.82 

3.74 

4.53 

.2 4.1 ± 1.2 2.55 ±0.25 Home. 

3 115 ±6 79 ±12 GALLEX 

) 115 ±6 74 ±16 SAGE 
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Table Illb: Characteristics of the DS Solar Models in Table Ilia 
^center; s=surface; bc=base of convective zone; N = log{[N]/[H] + 



Parameter 


DSND 


DS94 


DS95 


n [lO'^K] 


1.553 


1.554 


1.561 


Pc [g cmT^] 


154.9 


155.3 


155.4 




0.3491 


0.3462 


0.3424 


Yc 


0.6333 


0.6359 


0.6380 




0.01757 


0.01802 


0.01950 




0.6978 


0.7243 


0.7512 


Ys 


0.2850 


0.2597 


0.2308 


Zs 


0.01703 


0.01574 


0.0170 




8.55 


8.50 


8.55 




7.97 


7.92 


7.97 




8.87 


8.82 


8.87 




8.08 


8.03 


8.08 


X(>24 Mg)/Zs 


0.00464 


0.00414 


0.00415 


Rconv 


0.7306 


0.7105 


0.7301 


Tbc [lO^K] 


1.97 


2.10 


2.105 


Teff [K] 


5895 


5920 


5803 


^0 


0.69775 


0.69859 


0.7295 


Yo 


0.285007 


0.28267 


0.2509 


Zo 


0.01703 


0.01931 


0.01833 
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Table IVa: Comparison Between Solar Neutrino Fluxes Predicted 
by the Dar-Shaviv Models and the Bahcall-Pinsonneault best Solar 
Models. 



4>v{pp) [l^^^cm ^] 

M'^Be) [lO^cm-^s-i] 

M^B) [lO^cm-^s-i] 
M^^N) [10^cm-2s-i] 
M^^O) [lO^cm-^s-i] 

^{^a)ci [SNU] 
n^a)Ga [SNU] 



BP92 


DS94 


6.00 


6.06 


1.43 


1.42 


4.89 


4.00 


5.69 


2.60 


4.92 


3.30 


4.26 


3.19 


5.39 


3.84 


8± 1 


4.2 ± 1.2 


132 ±7 


116 ±6 



BP95 


DS95 


5.91 


6.10 


1.40 


1.43 


5.15 


3.71 


6.62 


2.49 


6.18 


3.82 


5.45 


3.74 


6.48 


4.53 



9.3 ±1.4 4.1 ±1.2 
137 ±8 115 ±6 
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Table IVb Characteristics of the BP95, DS94, and DS95 Solar 
Models in Table Ilia (c=center; s=surface; bc=base of convective 
zone; N = log([N]/[H]) + 12). 



Parameter 


BP95 


DS94 


DS95 


Tc [lO^K] 


1.584 


1.554 


1.561 


Pc [g cm~^] 


156.2 


155.3 


155.4 




0.3333 


0.3462 


0.3424 




0.6456 


0.6359 


0.6380 




0.0211 


0.01950 


0.01940 




0.7351 


0.7243 


0.7512 


Ys 


0.2470 


0.2597 


0.2308 




0.01798 


0.01574 


0.0170 




8.55 


8.50 


8.55 




7.97 


7.92 


7.97 




8.87 


8.82 


8.87 




8.08 


8.03 


8.08 


X(>24 Mg)/Z, 




0.00414 


0.00415 


Rconv 


0.712 


0.7105 


0.7301 


The [lO^K] 


2.20 


2.10 


2.105 


Teff [K] 




5920 


5803 
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Table V: Comparison between the solar neutrino fluxes calculated 
from the best standard solar model with no diffusion of BP95 and 
those calculated with the Dar-Shaviv SSM code with the same nuclear 
reaction rates, opacities, composition and astrophysical parameters. 
The predictions of the current best standard solar models of Dar & 
Shaviv and of Bahcall &: Pinsonneault are also included. 



BP95(ND) 



'u{pp) [W^^cm 


-2s-l] 


6.01 


'u{pep) [lO^cm" 




1.44 


'SCBe) [lO^cm 




4.53 




4.85 






4.07 


vC^O) [10^ cm" 




3.45 


i^(i^F) [lO^cm- 




4.02 



T.{(t>cT)ci[SNU] 7±1 
^{4>a)Ga[SNU] 127 ±6 



DS(BPND) BP95 DS95 



6.08 5.91 6.10 

1.43 1.40 1.43 

4.79 5.15 3.71 

5.07 6.62 2.49 

4.05 6.18 3.82 
3.38 5.45 3.74 

4.06 6.48 4.53 



7±1 9.3 ±1.4 4.1 ±1.2 

128 ±7 137 ±8 115 ±6 
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